Conductance of a helical edge liquid coupled to a magnetic impurity 
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Transport in an ideal two-dimensional quantum spin Hall device is dominated by the counterprop- 
agating edge states ol electrons with opposite spins, giving the universal value of the conductance, 
2e 2 fh. We study the effect on the conductance of a magnetic impurity, which can backscatter an 
electron from one edge state to the other. In the case of isotropic Kondo exchange we find that the 
correction to the electrical conductance caused by such an impurity vanishes in the dc limit, while 
the thermal conductance does acquire a finite correction due to the spin-flip backscattering. 

PACS numbers: 71.10.Pm, 72.10.Fk 



CO 



T3 
C 

o 

> 
CN 

CN 
O 



X 



Topological insulators have been actively studied in the 
last few years [l[ . In these systems the bulk of the sample 
is insulating, with gapless electronic excitations allowed 
only at the boundary. One example of this phenomenon, 
the quantum Hall effect, is realized in two-dimensional 
electron systems placed in a strong magnetic field per- 
pendicular to the sample. Another example is the so- 
called quantum spin Hall (QSH) effect, recently observed 
in HgTc quantum wells [2[ . It is present in time-reversal 
invariant two-dimensional systems, and can be viewed as 
two coexisting quantum Hall states for opposite spin sub- 
systems. Although no external magnetic field is present, 
an effective field is generated via spin-orbit coupling. The 
sign of the field is opposite for spin-f and spin-| electrons, 
resulting in two edge states propagating in opposite di- 
rections. 

In the absence of defects, the conductance measured 
between two contacts attached to the edges of a quantum 
spin Hall device is 2e 2 /h. (Here e is the electron charge 
and h is the Planck's constant.) This can be understood, 
e.g., as twice the conductance of a standard quantum Hall 
structure, as the number of the edge states is doubled. 
The presence of two countcrpropagating edge states also 
allows for the possibility of backscattering of electrons, 
which would reduce the conductance. However, such 
backscattering processes are strongly restricted by the 
time-reversal symmetry 0,0]. In particular, an impurity 
without internal degrees of freedom cannot backscatter a 
single electron at the edge, though backscattering of two 
electrons is allowed. Another allowed process is backscat- 
tering of a single electron by a magnetic impurity. 

The effect of a magnetic impurity on the conductance 
of a QSH device was discussed recently by Macicjko et al. 
[j| . They argued that at high temperature the backscat- 
tering by magnetic impurity gives a small negative cor- 
rection to the conductance. As the temperature T is low- 
ered, the Kondo effect enhances the backscattering and 
the correction grows. The trend then reverses at T of the 
order of the Kondo temperature, and at T = the quan- 
tized conductance 2e 2 /h is restored. In this Letter we 
study the correction SG to the conductance as a function 



of frequency uj, focusing on the case of isotropic Kondo 
exchange and neglecting the two-particle backscattering. 
Our results for 6G(u>) agree with the predictions by Ma- 
ciejko et al. [B| when ui exceeds a certain relaxation rate, 
whereas in the dc limit u) — > we find 6G = 0. 

Following Refs. we combine the two edge states 

into a single nonchiral Tomonaga-Luttinger liquid and 
write its Hamiltonian in the standard bosonized form IfJ 



H 



hv 
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Here v is the velocity of the edge states, and parameter 
K accounts for the interactions between the electrons at 
the edge, with the repulsive interactions corresponding 
to K < 1. The bosonic fields 6(x) and <f>(x) satisfy the 
commutation relations [4>(x),d y 0(y)} = m8(x — y). For 
simplicity, we limit ourselves to the case of spin-i im- 
purity, and write the Kondo coupling to the edge states 
using the standard boson representation of the electron 
operators f/'t.i = (2na)~ 1 ^ 2 e~^ e± '^\ where a is the short 
distance cutoff and we assigned spins up and down to the 
right- and left-moving electrons, respectively We allow 
for the possibility of anisotropic coupling and write the 
spin-flip term as 
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-2i<j>(0) g- 
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where S = S x ± iS y and S l are spin operators. In 
addition, the z component of the impurity spin is coupled 
to the Tomonaga-Luttinger liquid via 



H- 



-d x 9(0)S z 

7T 



(3) 



For convenience, our Hamiltonian allows for uniaxial ex- 
change anisotropy, which disappears at J z = J±. Up to 
a minor difference in notation, the Hamiltonian {l])-(|3]) 
coincides with the expression used in Refs. [H, H[- 

To discuss the conductance of the device, one has to 
add to the Hamiltonian a term accounting for the applied 
bias. The most straightforward approach is to assign 
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different chemicai potentials ±eV/2 to the right- and left- 
moving electrons. Such a perturbation takes the form 



eV 

H v = -— I dxd x 9. 



(4) 



One should note that assigning the same chemical poten- 
tial to all the electrons moving in the same direction is 
an approximation. It assumes that the backscattering by 
impurity is weak, SG <C e 2 /h, and that ui <C v/L, where 
L is the distance between the source and drain contacts. 

The operator (01 commutes with Hq and H z , and its 
only effect is to change the energy by ±eV each time an 
electron is backscattered and the impurity spin is flipped 
by H±. The same effect is achieved by assigning dif- 
ferent energy values to the up and down components of 
the impurity spin, i.e., by replacing (@| with an effective 
magnetic field term 



H v = -eVS z . 



(5) 



This replacement can be more formally justified by notic- 
ing that the difference of the operators (j4|) and ([5]) com- 
mutes with the Hamiltonian. 

To simplify the subsequent calculations it is convenient 
to rescale the bosonic fields {(f>, 6} — s- {\/K<p, 6/y/K} and 
then perform the unitary transformation U = e lX( t>(°) s * . 
For A = Jz/ttHv^/K the additional term arising from the 
transformation of Hq cancels the coupling ([3]) of the z 
components of spins, resulting in the Hamiltonian 



Ho = ^ I dx[{d x 6f + {d^f], 



(6) 



H ± = ^[ e ^0(o)s+ + e -*2^(o) 5 - ]; (7) 
2na 

with the bias contribution retaining its form. The 
advantage of this procedure is that the effect of electron- 
electron interactions and coupling of z components of 
spins are now accounted for by a single parameter K = 
K{1 - J z /2irhvK) 2 . 

Our first goal is to evaluate the correction 8G(us) to 
the conductance of the system due to the backscatter- 
ing of electrons by the impurity. Since backscattering 
of a right-moving electron is always accompanied by the 
impurity spin flip from J, to f , the correction to the cur- 
rent can be computed as the time derivative of S z , i.e., 
SI = -ed t S z . With the Hamiltonian in the form ©-© 
one immediately obtains 



SI= — Ij_\ e ^V&4>(0) g+ 
h 2ira l 
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In the linear response theory the conductance can be 
found using the Kubo formula, which expresses SG(u>) 
in terms of the current-current correlator. The latter 
cannot be found exactly for arbitrary values of the pa- 
rameters K and J±. Assuming weak coupling ([7]). one 



can evaluate SG to the lowest order in J±, which results 
in the following expression: 

2e 2 f J\\ 2 fnT\ 2k 



(9) 



SG ™ = -W{2^ (-D Sin( ^ 
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Here we have introduced the bandwidth D = hv/a. 

In the most interesting regime of low frequencies, fuv <C 
T, the conductance is 



SG = 



e 2 7o 
2T ' 



where 70, defined as 



7o = JlT, 



[T(K)} 2 /2^T\ 
(2nh) 2 avT(2K) \ D ) 



2K-1 



(10) 



(11) 



has the meaning of the rate of spin-flip processes at V = 
0. The same result was obtained by a similar method 
by Maciejko et al. It is important to realize that the 
derivation relied on the perturbation theory in J±. Let 
us now show that at low frequencies w < 70 « J| the 
conductance deviates from (fT0|) . 

To this end we consider the dynamics of the impurity 
spin at finite temperature and weak spin-flip coupling, 
70 <C T/h. In this case, the spin remains in the j" or I 
state for a long time. The rates of spin-flip events are 
easily computed and take the form 



7± = 7o 



2T 



(12) 



for small voltage eV -C T, where 7-1- and 7_ are the rates 
of up-flip and down-flip, respectively. 

At any given moment the impurity spin can be in ei- 
ther j" or I state, and its behavior is described by the 
probabilities P-j- and P±. Their time dependence can be 
found from a simple rate equation 



d t P\ = 7+P4 - 7_Pf 



(13) 



and the condition Pj- + Pj_ = 1. Each spin flip is accompa- 
nied by backscattering of a single electron. The resulting 



correction to the electric current is 81 



edtPf. To 



find the linear conductance one substitutes the rates in 
the form (12]) with V = V e~ lujt into ([13]). and calculates 
SG = 8I/{V G e~ iu}t ). This procedure yields 



SG(lu) = 



e 2 7o 



2T oj + 2i-f 



(14) 



The rate equation approach is applicable only at rela- 
tively low frequencies, fku <C T. In a broad frequency 
range 70 <C w <C T/h expression p^|) recovers the per- 
turbative result ^U|) . 
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The results (jlOl) and (|14l) differ dramatically in the dc 
limit, oj — y 0, where the correction (|14[) vanishes. This is 
easily understood if one notices that every time an elec- 
tron is backscattered to the left (right) the impurity spin 
is flipped up (down). Thus, the backscattering current 
changes its direction with every spin flip, and in the dc 
limit the correction to the conductance vanishes. This ar- 
gument also applies at T < fvyo, when the rate equation 
approach is no longer applicable. 

To illustrate this point, we consider a special case K = 
i, in which the dimension of the operator exp[i2\//F<^(0)] 
in the definition of the spin- flip operator ([7} is |, [5j. 
This case corresponds to the well-understood Toulouse 
limit of the Kondo model, in which the Hamiltonian can 
be reduced to that of a system of noninteracting chiral 
spinless fermions ^(x) coupled to a discrete level: 



H = -ihv I ¥{x)d x ^{x)dx-eVSd 
J± 



= [d t *(0) + * t (0)d]. 



(15) 



Here d is the annihilation operator of a fermion at the 
discrete level which models the impurity spin, S z = d* d— 
■i. The model (fl~5j) is easily solved exactly, resulting in 



5G(u) 



hio-^+b+iS (£! 2 +r 2 )(£ 2 2 +r2) 



(16) 

Here the level width T — J^/Airhav coincides with hjo 
for K = 1/2, Eq. {HJ, and = l/(e« /T + 1) is the 
occupation probability of a fermion state with energy £ 
measured from the Fermi level. 

In the high-temperature limit T>T one can expand 
the difference of the Fermi functions in the numerator 
to linear order in £2 — £i- This approximation recovers 
our earlier result (fT4|) with 70 = T/H. Importantly, the 
Toulouse limit expression (fTtJj) for 5G(w) does not rely 
on the rate equation approach and is therefore valid for 
any temperature. It is easy to see that in the dc limit 
the expression (fT6|) vanishes for any T. The dependence 
SG(u>) at various temperatures is shown in Fig. [T] 

In the absence of corrections to the quantized dc 
conductance of the system, it is interesting to explore 
whether other transport properties are affected by cou- 
pling to the impurity. One such property is the ther- 
mal conductance JC. Let us assume that the left- and 
right-moving electrons originate in the leads with differ- 
ent temperatures, T and T + AT, respectively. Spin-flip 
scattering by the impurity gives rise to two kinds of elec- 
tron scattering processes. The electrons are either scat- 
tered from the warmer right-moving system to a colder 
left-moving one, or in the opposite direction. These two 
types of processes change the charge of the right-moving 
system by e and — e, respectively, resulting in no correc- 
tion to the time-averaged electric current. On the other 
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FIG. 1: (color online) Real part and imaginary part of 8G(uj) 
for several temperatures T/T. 



hand, each scattering process transfers heat in the same 
direction: from the warm right-moving system to the cold 
left-moving one. As a result, one expects to find a finite 
negative correction SfC to the thermal conductance of the 
system. 

To evaluate SfC, we identify the operator of the differ- 
ence of energies in the right- and left-moving branches 
SH = (Hff - Hq)/2. The energy densities of the two 



subsystems are given by Hv[d x 



hv 



1 T 1 



' /Aix, leading to 



SH = —— I d x (j)d x 9dx. 
2n 



(17) 



The operator SJe of the energy current transferred be- 
tween the right- and left-moving branches is then found 
as the time derivative of p7|) . resulting in 



SJ E 



■d x (t>(0)SI, 



(18) 



where SI is the electric current operator in Eq. ©. 

It is worth noting that the operators d x (j>(0) and 51 act 
in separate subspaces of odd and even <p(x). In particular, 
the dynamics of d x (j)(0) is not affected by the coupling 
([7]) to the impurity. A similar observation was made by 
Kane and Fisher [7|] in the study of thermal transport of a 
Luttinger liquid through a tunneling barrier. They used 
the factorization of the energy current operator similar 
to Eq. (fl~8|) to obtain a relation between the thermal and 
electrical conductances of the system. Repeating their 
procedure for our Hamiltonian, we obtain the expression 



SfC 



h 3 K 



oj 2 RcSG(u}) 



8e 2 T 2 J_ oc sinh 2 (nw/2T) 



du, 



(19) 



fully analogous to Eq. (19) of Ref. Q. 

The relation (p~9|) shows that thermal conductance is 
determined by the electrical conductance at frequencies 
u> ~ T/H. As a result, the fact that SG vanishes at u — > 
does not mean that there will be no correction to JC. For 
instance, in the lowest order in spin- flip scattering (J7J, 
one can find SfC by substituting into Eq. (|T9"|) SG{uS) in 
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FIG. 2: Nonmonotonic temperature dependence of the ther- 
mal conductance K, — ICo + 5IC in the Toulouse limit, Eq. (|21|l . 



the form ([9]). This procedure yields 



5K = -- 



K 2 



2 2K+1 



To- 



(20) 



Within the applicability of the perturbation theory, 
tvjo *C T, this correction is small compared to the nomi- 
nal value /Co = tt 2 T/3Ii of the thermal conductance of a 
Kramers' pair of edge states. 

At ~ T the perturbative approach is not appli- 
cable, but an exact solution is possible in the Toulouse 
limit, K = 1/2. Substituting the conductance (fT6|) into 
Eq. (fT9|) we obtain 



SK = 



fe-6) 3 K6)-^(6)]^irf6 

HhT* J J (g + + I*) sinh 2 [(6 - 6)/2T] ' 

(21) 

At high temperatures (T > T), <5/C = -7r 2 r/16?i, in 
agreement with the perturbative result (|2"0")) . On the 
other hand, the fact that Re5G(u) oc u 2 at w — > means 
that the correction <WC in Eq. (TIT)]) is suppressed as T 3 at 
low temperatures, T <C T. Indeed, from Eq. (|2"Tj) one 
finds SIC = -2TT 3 T 3 /l5hT 2 in this regime. 

The suppression of SIC in both the low- and high- 
temperature limits results in nonmonotonic behavior of 
the normalized thermal conductance IC/ICq as a function 
of temperature, Fig.[2J Such nonmonotonic behavior was 
predicted for the electrical conductance by Maciejko et 
al. Q , based on the well-known nonmonotonic tempera- 
ture dependence of the spin-flip scattering in the Kondo 
problem. Although our theory predicts quantized dc con- 
ductance 2e 2 /h at any temperature, the nonmonotonic 
temperature dependence is recovered for thermal trans- 
port. 

At T < 70 and uj ~ 70 our correction to the electrical 
conductance can take rather large values |(5G(w)| ~ e 2 /h, 
see, e.g., Fig. [TJ Wc do not believe the correction 5G(uj) 
has a clear physical meaning in this regime, as the condi- 
tions for the applicability of the model (j4j are violated. 
This caveat, however, does not apply to our main con- 
clusion, namely vanishing SG at uj — > 0. Similarly, the 



behavior of the thermal conductance shown in Fig. [2] is 
only qualitatively correct at T ~ T, where the relative 
correction SIC/ICq is of order unity. 

The main feature of our work, the absence of correc- 
tion to the dc conductance due to Kondo scattering, is 
caused by the symmetry of the Hamiltonian ([IJ-Q that 
preserves the z component of the total spin of the sys- 
tem. Our model is justified if the magnetic impurity is 
approximated as a single-orbital Anderson impurity, in 
which case the Kondo exchange is isotropic, Jcr ■ S, even 
in the presence of spin-orbit coupling |8|. Anisotropic 
corrections may appear in the case of multiple orbitals 
Exchange anisotropy can break the conservation of 
the z component of the total spin and result in nonvan- 
ishing correction to the dc conductance. In the simplest 
model of anisotropic exchange Ji^S 1 with J x ^ J y the 
correction to the dc conductance above the Kondo tem- 
perature is easily found from the rate equation approach, 



2 1 2 



SG 



AT 31 



J 2 



(22) 



with T defined by Eq. (fTTjl. At T below the Kondo tem- 
perature the impurity spin is screened, and the effect of 
exchange anisotropy is reduced to that of an impurity- 
induced two-electron backscattering. The effect of these 
processes on conductance was studied in Ref. |f| and 
found to be suppressed at T -> as T 2 ! 4 ^- 1 ). 

Another limitation of our model is the assumption that 
all electrons moving in the same direction have the same 
spin. Kramers degeneracy guarantees that electrons with 
momenta ±p have opposite spins, and wc defined the 
z axis as the direction of the spin of the right mover 
at the Fermi level, which is determined by the specific 
form of spin-orbit coupling. Transport is controlled by 
electrons with momenta |p±pi?| ~ T/v, whose spins 
may deviate slightly from the z direction defined at p = 
pp. This deviation may result in a correction to the dc 
conductance vanishing at low temperatures as a power of 
T/D. 
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